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Abstract 

We prove the global existence of weak solutions for the 2-D compressible Navier-Stokes equations 
with a density-dependent viscosity coefficient (A = A(p)). Initial data and solutions are small in 
energy-norm with nonnegative densities having arbitrarily large sup-norm. Then, we show that if 
there is a vacuum domain at the initial time, then the vacuum domain will retain for all time, and 
vanishes as time goes to infinity. At last, we show that the condition of p =constant will induce 
a singularity of the system at vacuum. Thus, the viscosity coefficient p plays a key role in the 
Navier-Stokes equations. 

1 Introduction 

In this paper, we consider the following 2-D compressible Navier-Stokes equations 

p t + div(pu) = 0, . , 

(pu) t + div(pu (g) u) + VP = ^Am + V((/U + A(p))divu) + pf, 1 ' ' 

for i£l 2 and t > 0, with the boundary and initial conditions 

u(x,t) — > 0, p(x,t) — > p > 0, as \x\ — > ao, t > 0, (1.2) 

(p,u)\ t= o = (p ,u ). (1.3) 

Here p(x,t), u(x,t) and P — P{p) stand for the fluid density, velocity and pressure respectively, / is 
a given external force, the dynamic viscosity coefficient p is a positive constant, the second viscosity 
coefficient A = X(p) is a function of p. 

At first, we prove the global existence of weak solutions that are in an "intermediate" regularity class 
in which energies are small, but oscillations are arbitrarily large. Specifically, we fix a positive constant p, 
assume that (po — P, uq) are small in L 2 , and (po, uq) £ L°° x H 1 with no restrictions on their norms. Our 
existence result accommodates a wide class of pressures P, including pressures that are not monotone 
in p. Since the solutions may exhibit vacuum states and discontinuities in density and velocity gradient 
across hypersurfaces, our results are consequently much less regular and much more general than the 
well-known small-smooth theory, such as [H Q~5] . This existence result generalizes and improves upon 
the earlier result of Vaigant-Kazhikhov [16] in two significant ways: the space domain is unbounded and 
the initial density may vanish in an open set. It also generalizes and improves upon earlier results of 
Hoff-Santos [6] and Hoff [3 [9] in two significant ways: the second viscosity coefficient A is a function 
of the density p, and we omit the condition J(l + |x| 2 ) a (po| w o| 2 + G(po))dx < Mo with a constant a > 0. 

We now give a precise formulation of our existence result. 

Definition 1.1. We say that (p,u) is a weak solution of (|l.ip (|1.3[) . if p and u are suitably integrable 
and satisfy that 
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p<pdx 



for all times t 2 > h > and all 4> G Cq(K 2 x [ti,t 2 ]), 



ti Jii 
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64 + pu ■ \/<j))dxdt 



(1.4) 



*2 /-t 2 



{pit • "0* + p(it • V^) • it + Pdiv-0 }dxdt 



ti -'ti 

{pdkU 3 dkip 3 + (p + A)divudiv^ — pfip}dxdt 



puipdx 

'ii 

for all times t 2 > *i > and all ip = (i/j 1 ,^ 2 ) € (Cq (M 2 x i 2 ])) 2 - 
Concerning the pressure P, viscosity coefficients p and A, we fix < p < p and assume that 

( PeC\[0 : p}), AeC 2 ([0,p]), 
P>0, A(p)>0, pe [0,p], 
P(0) = 0, P'(p) > 0, 

(p - p)[P(p) - P(p)] >0,^ft P e[0, p], 

PeC 2 ([0,p\) or 2^A(.) is a monotone function on [0, p]. 
Let G be the potential energy density, defined by 

G{p) = p J p (' s ) _2LL ds. 



(1.5) 



(1.6) 



(1.7) 



Then for any g 6 C 2 ([0,p]) with g(p) = g'(p) = 0, there is a constant C such that \g(p)\ < CG(p), 
pe [0,p]. 

We measure the sizes of the initial data and the external force by 



Co = 



Po\u \ 2 + G(p ) dx 



C / =sup||/(-,t)||| 2 + / (||/(,i)|| i2 + ||/(.,<)|| 2 2 +a 4 ||V/||i 4 +a 3 ||/,(,t)||| 2 )A, 



(1.8) 
(1.9) 



and 



M q = ^(^ 2 ||/ t ||i 2 +^ 3 ||V/||t 4 +^ + *||/ t ||^ e +a 2 +4||V/||^) 

+ / |Vw | 2 da; + sup||/|| L2+g , 
J t>o 

where a(t) = min{l,t} and q is a constant satisfying 



dt 



(1.10) 



<7 G (0, 2) and q 2 < 



4p 



Vpe [0,p]. 



P + A(p) 

As in [SI [7] , we recall the definition of the vorticity matrix w 1,k = d^v? —djU k , and define the function 

F= (A + 2p)divu-P(p) + P(p). (1.11) 

We also define the convective derivative ^ by ^ = li; = 1^ + u • Vw, the Holder norm 

\v(x) — u(ty)| 
<v>%= sup J-U 
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and 

. s c,i9 _ \g(x,t) -g{y,s)\ 

^ 9 >Ax[t x ,t a ]— . ,, b N up „ r . U_,,|at -L |/_ ,1/3' 

(a;,t)/(y,s) 

where o : i C I 2 ^ R 2 , g : A x [t u t 2 ] -> R 2 and a,/3 G (0, 1]. 
The following is the global existence result of this paper. 

Theorem 1.1. Assume that conditions il.6\) - !!l.l0\) hold. Then, for a given positive number M (not 
necessarily small) and pi £ (p, p), there are positive numbers e, C and 9, such that, the Cauchy problem 
11. 1) ) - fO)) with the initial data (po,wo) and external force f satisfying 

< inf p < sup p < pi , 

C + C f <e, (1.12) 
M q < M, 



has a global weak solution (p, w) in the sense of |i.^| )-( T7T5p satisfying 

C _1 inf po < P < P, a.e. (1.13) 

j p-p, pu g cao,^)-^- 1 ), . . 

\ Vn€i 2 (l 2 x [0,00)), 1 j 

^^oo)+ su P(l|V^(^)IU 2 + ||V W (.,t)|| L2 )<C(a,r)(C + C' / ) fl , (1.15) 



t>T 



sup / |Vu| 2 dx + / / a\Vu\ 2 dxds < C, (1.16) 
t>0 J Jo J 

T 

\\F(-,t)\\ L -dt<C(T), (1.17) 



< ^ >R»^?1 + < w >R ™^ C ^ T )' (1-18) 
where a G (0, 1), r > and a' G (0, arl ^ 

sup y QpI^I 2 + IP - Pi 2 + cr|Vw| 2 ^ dx 

+ y y (|Vw| 2 + a\(pu) t + div(pw (& u)\ 2 + cr 2 |V?i| 2 ) dxdt 
< C(C + Cf) e . (1.19) 
In addition, in the case that inf po > 0, the term J °° f a \u\ 2 dxdt may be included on the left hand side of 

Remark 1.1. Here, 9 is a universal positive constant (we choose 9 = i in this paper), e and C depend on 
p, pi, p, P, A, p, q and M. 

Remark 1.2. For example, we can choose that P = ^4p 7 and A(p) = cp' 3 with 7 > 1 and /3 > 2, where A 
and c are two positive constants. Also, we can choose that A is a positive constant and P = Ap 1 with 
7 > 1. 

Remark 1.3. Considering the non- vacuum case, i.e., po > 2p > 0, we can replace the condition (|1.6|) by 



f PeC^O.p]), AeC%p]), 
p>0, A( P )>0, P e[0,p] 

P(0) = 0, P'(p) > 0, (1.20) 
(p - p) [P(p) - P(p)] > 0, p ^ p, p G [p, p] , 
P G C' 2 ([p, p]) or 9^^.) is a monotone function on [p, p], 

where < 2p < p < p. Then, we can choose that P = Ap 1 and A(p) = cp' 3 with 7 > 1 and (3 > 0. 
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Remark 1.4. Using a similar argument as that in [16], one can obtain the uniqueness of the solution if 
the initial data satisfy po > p > 0, and (po, u ) £ C 1+q x C 2+a or (p , u ) £ W x - q x H 2 , where a £ (0, 1) 
and q > 2. 

The proof of Theorem 11.11 consists in the derivation of a priori estimates for smooth solutions cor- 
responding to mollified initial data, and the application of these estimates in extracting limiting weak 
solutions as the mollifying parameter goes to zero. Specifically, in Section [21 we fix a smooth, local in 
time solution for which < p < p and A x + A 2 < 2(C + C f ) e ', where 9 £ (0, 1), 



and 



Ai{T) = sup a [\\7u\ 2 dx + [ [ ap\u\ 2 dxdt 

0<t<T J JO J 

A 2 (T) = sup a 2 [ p\u\ 2 dx + [ [ a 2 \X7u\ 2 dxdt, 

0<t<T J JO J 



then obtain the estimate A\ + A 2 < (Co + C/) e , and prove that the density remains in a compact subset 
of [0,/5). Using the classical continuation method, we can close these estimates, which are then applied 
in Section [3] to show the solution can be obtained in the limit as the mollifying parameter goes to zero. 

Using the initial condition uq £ H\ we can obtain pointwise bounds for F in Proposition 12.61 which 
is the key point of the a priori estimates. Because that the mass equation can be transformed to the 
following form, 

j t A(p(x(t),t)) + P( P (x(t), t)) - P(p) = -F(x(t),t), (1.21) 

where A satisfies that A(p) = and A'(p) — HM^M^ a curve x(t) satisfies x(t) — u(x(t), t), thus pointwise 
bounds for the density will therefore follow from pointwise bounds for F. On the other hand, using a 
similar argument as that in [16] and the estimate J„ \\F(-,t)\\Lx>dt < C(T), we can obtain the strong 
limit of approximate densities {//}, see Section [3] 

Then, we study the propagation of singularities in solutions obtained in Theorem 11.11 We show that 
each point of M 2 determines a unique integral curve of the velocity field at the initial time t = 0, and that 
this system of integral curves defines a locally bi-H61der homeomorphism of any open subset f2 onto its 
image f2* at each time t > 0. Using this Lagrangean structure, we show that if there is a vacuum domain 
at the initial time, then the vacuum domain will exist for all time, and vanishes as time goes to infinity, 
see Theorem 11.41 Also, we show that, if the initial density has a limit at a point from a given side of 
a continuous hypersurface, then at each later time both the density and the divergence of the velocity 
have limits at the transported point from the corresponding side of the transported hypersurface, which 
is also a continuous manifold. If the limits from both sides exist, then the Rankine-Hugoniot conditions 
hold in a strict pointwise sense, showing that the jump in the (A + 2p)divu is proportional to the jump in 
the pressure. This leads to a derivation of an explicit representation for the strength of the jump in A(p) 
in non-vacuum domain. These results generalize and improve upon the earlier results of Hoff-Santos [6] 
in a significant way: the domain may contain the vacuum states. 



Theorem 1.2. Assume that the conditions of Theorem \l.l\ hold. 

(1) For each to > and each xo £M. 2 , there is a unique curve X(-;xo,to) £ C 1 ((0, oo); R 2 )nC 1 ~'^ ([0, oo); R 2 ), 

a £ (0, 1), satisfying 

X(t; x , t ) = x Q + / u(X(s;x ,t ),s)ds. (1-22) 
J to 

(2) Denote X(t,Xo) — X(t;xo,0). For each t > and any open set fl C R 2 , 11* = X(t, -)f2 is open and 

the map xo ' — ► X(t,xo) is a homeomorphism of O onto f2'. 

(3) For any < ti,t2 < T, the map X(t\,y) — > X{t 2 ,y) is Holder continuous from R 2 onto R 2 . 

Specifically, for any y\,y 2 £ R 2 , 

\X(t 2 , y 2 ) - X(t 2 , yi )\ < exp(l - e- c ^)\X{tx,y 2 ) - X^y^^ . (1.23) 
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(4) Let M C R 2 be a C a 1-manifold, where a G [0, 1). Then, for any t > 0, M 1 = X(t, -)M is a C 13 
1-manifold, where [3 = ae~ c ^ 1+t \ 



Theorem 1.3. Assume that the conditions of Theorem \l.l\ hold. Let V be a nonempty open set in R 2 . 
// essinfpol v > P > 0, then there is a positive number p~ such that, 

p(-,t)\v< > £", 

for all t > 0, where V 1 = X(t, -)V . 

Theorem 1.4. Assume that the conditions of Theorem \l.l\ hold. Let U be a nonempty open set in R 2 . 
Assume that p 1 f/ = 0. Then, 

p{-,t)\ ut = 0, 

for all t > 0, where U = X(t, -)U . Furthermore, we have 

lim / \p- p\ A {x,t)dx = 0, (1.24) 



t — >oo 

and 

lim \{x e R 2 \p(x, t) = 0}| = 0. (1.25) 

t — >oo 

Theorems 11.2111.41 are proved in Section [4j 

In the following theorem, applying the Lagrangean structure of Theorem 11.21 we establish a result 
concerning the transport by the velocity field of pointwise continuity of the density. Recall first that the 
oscillation of g at x with respect to E is defined by (as in [5]) 

osc(s; x, E) = lim {esss\xpg\ E nB R {x) ~ ess\nig\ ErtBR{x) ) , 
ii — >o 

where x S E and g maps an open set £cl 2 into R. We shall say that g is continuous at an interior 
point x of E, if osc(g; x, E) = 0. 

Theorem 1.5. Assume that the conditions of Theorem li.il hold. Let E C R 2 be open and xq G E. If 
osc(po; xq, E) — 0, then osc(p(-, f); X(t, Xq), X(t, -)E) = 0. In particular, if Xq S E and p$ is continuous 
at xq, then p(-,t) is continuous at X(t,Xo). 

Now, let M. be a C° 1-manifold in R 2 and xq G M.. Then there is a neighborhood G of xq which 
is the disjoint union G = (G n M) U ^+ U where £7± are open and Xo is a limit point of each. If 
osc(g; xo, E + ) — 0, then the common value g{xo+, t) is the one-sided limit of g at xq from the plus-side of 
M., and similar for the one-sided limit g(xQ—,t) from the minus-side of M.. If both of these limits exist, 
then the difference [g(xo)] := 9i x o+) — 9(xq—) is the jump in g at x n with respect to M. (see [5]) 

Now, we state our main results on the propagation of singularities in solutions. 

Theorem 1.6. Let (p,u) as in Theorem \l.l\ M. be a C° 1-manifold and xq G M.. 

(a) If pa has a one-sided limit at xq from the plus- side of M., then for eacht > 0, p(-,t) and divu(-,i) 
have one-sided limits at X(t,Xo) from the plus-side of the C° 1-manifold X(t,-)A4 corresponding to the 
choice E\ = X(t,-)E + . The map t i— * p(X(t,xo)+,t) is in ([0, oo)) n C 1 ((0,oo)) and the map 
t i ► drvu(X(t, x ) + , t) is locally Holder continuous on (0,oo). 

(b) If both one-sided limits pq(xq±) of po at Xq with respect to M exist, then for each t > 0, the jumps 
in P(p(-,t)) and divu(-,i) at X(t,Xo) satisfy the Rankine-Hugonoit condition 

{(2p + X(p(X(t, x Q ),t)))divu(X(t, x ), t)] = [P(p(X(t, x ), t))]. (1.26) 

(c) Furthermore, if po(xo±) > p > 0, then the jump in A(p) satisfies the representation 

[A(p(X(t,x ),t))} = exp (- f a(T,x Q )dT) [A(p (x ))} (1.27) 



where a(t xn) = [p(p(x(t ' ao) ' t)}1 
wriere a{i,,x ) - [A( p (x(t,x ),t))] 
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Theorems I1.5H1.6I are proved in Section 

Remark 1.5. From Theorems I1.2IILB1 we have that if [p ] is nonzero at Xq, then [/>(■, t)] is nonzero at 
X(t; xo, 0) for every t > 0. That is, singularities of this type persist for all time. On the other hand, if 
P'{p) > with p~ < p < ~p, then a is strictly positive and the jump in A(p) in the non-vacuum domain 
decays exponentially in time. Thus, the jumps in p, in x(p)+2p an< ^ m divu i n the non-vacuum domain 
decay exponentially in time as well. 

At last, we will show that the condition of p =constant will induce a singularity of the system at 
vacuum in the following two aspects: 1) considering the special case where two fluid regions initially 
separated by a vacuum region, we show that the solution we obtained is a nonphysical weak solution in 
which separate kinetic energies of the two fluids need not to be conserved; 2) we show the blow-up of 
smooth solutions for the spherically symmetric system when the initial density is compactly supported. 
Thus, the viscosity coefficient p plays a key role in the Navier-Stokes equations. 

If we consider the special case where two fluid regions initially separated by a vacuum region, Theorem 
11.41 shows that the solution obtained in Theorem 11.11 is a nonphysical weak solution in which the two 
fluids cannot collide independent of their initial velocities. In the following, we will show that the separate 
kinetic energies of the two fluids needn't to be conserved. 

If the initial data are spherically symmetric, i.e, 

x 

p (x) = g (r), u (x) = voir)-, r = \x\, (1.28) 

r 

and the external force / = 0, from Theorem ll.il one can prove that the system (|l.ip - (11.3j) has a spherically 
symmetric solution (p, u) satisfying 

p{x,t) = g{r,t), u( X> t)=v(r,t)^, r = \x\. (1.29) 



Then (g, v) is a solution of the following system 



d t g + d r (gv) + ±gv = 0, 

g(d t v + vd r v) + d r P{g) = d r {(X(g) + 2p)(8 r v + f )}, 



(1.30) 



v(r, t) -> 0, g(r, t) -»• p > 0, as r -> oo, t > 0, (1.31) 
(g,v)\ t= o = {go,v ). (1.32) 
Furthermore, assume that there are two positive constant < a < b, such that 

g (r) = 0, r e (a,b), 

and 

Qo(r) > g > 0, re (0, a) U (b, +oo). 
Then, from Theorems 1 1 . 2lfl~4l we have there are two curves a(t) and b(t) satisfying 

< a(t) < b(t) < oo, 

a(0) = a, a(t) = v(a(t), t), 6(0) = b, b'(t) = v(b(t), t), 
g(r,t) = 0, re (a(t),b(t)), 

and 

g(r, t)>g~>0, re (0, a(t)) U (b(t), +oo), 

for some positive constant g~ . Using a similar argument as that in [10j . we can obtain the following 
theorem. 



G 



Theorem 1.7. Assume that the conditions of Theorem \1.1\ hold, s 2 P(s)ds < oo, / = 0, and the 
initial data satisfy \1.28\) , then we have 

^S(t) = 2(A(0) + 2fj,)a{t)v{a{t),t) fl2( -^ _ &2 ( t ) ' f - °> ( L33 ) 



/•<»(*) /x _ \ r t r a ( s ) v 

E(t)=J {-gv 2 + G(g))rdr + J J (A + 2p)(v r + -) 2 rdrds, 

and 

— f p P(s) 

G(p)=P A^ds. 
Jo s 

Remark 1.6. If the viscosity coefficient p is a function of the density and A(0) = p(0) = 0, the equality 
(| 1 . 33[) implies the separate kinetic energies of the two fluids are conserved. Thus, the main reason for 
the appearance of non-physical solutions comes from the viscosity coefficient p being independent of the 
density. 

Remark 1.7. The physical solution of the (|1.30p - (|1.32[) may be obtained by constructing separately the 
solutions for each of the fluids £>o|[o,a] an d Qo\\b,oo) with the boundary conditions 

(A(f?) + 2p)(v r + -)= P(g), r = a{t), b(t). (1.34) 
r 

When a(t) < b(t), one can obtain the composite solution (g,v). Observe that the kinetic energies are 
separately conserved because of the boundary conditions (|1.34[) . When vq is large and positive on [0,a], 
large and negative on [6, oo), a collision a(t) — b(t) may occur in finite time. 

Finally, in Section we will give a non-global existence theorem on smooth solutions for the spher- 
ically symmetric system when the initial density is of compact support. The corresponding theorem 
on compressible Navier-Stokes equations with constant viscosity and heat conductivity coefficients was 
obtained in [17] • Here we generalize the above theorem to the case when the second viscosity coefficient 
depends on the density for the isentropic gas flow. 

Theorem 1.8. Suppose that (p,u) € C 1 ([0,T\;H k ), k > 3 is a spherically symmetric solution to the 
Cauchy problem U.l\) and U.3\) with f — 0. Assume that P(p) = Ap 1 and A(p) = cp@ with 1 < [3 < 7 
and A, c > 0. If the support of the initial density po is compact and po ^ 0, then T must be finite. 

We now briefly review some previous works about the Navier-Stokes equations with density-dependent 
viscosity coefficients. For the free boundary problem of one-dimensional or spherically symmetric isen- 
tropic fluids, there are many works, please see [SJ [HI HH [TU [19] and the references cited therein. Under 
a special condition between p and A, A = 2/3/i' — 2/z, there are some existence results of global weak 
solutions for the system with the Korteweg stress tensor or the additional quadratic friction term, see 
[lj [2] ■ H. L. Li, J. Li and Z. P. Xin [12] showed a very interesting result that for any global entropy weak 
solution of the one-dimensional system, any vacuum state must vanish within finite time. Also see Lions 
[13] for multidimensional isentropic fluids. 

We should mention that the methods introduced by Hoff-Santos in [5], Hoff in [S] and Vaigant- 
Kazhikhov in [16] will play a crucial role in our proof here. 



2 A priori estimates 

In this section, we derive some a priori estimates for local smooth solutions of the system (|l.l[) - (|1.3p with 
strictly positive densities. Thus, we fix a smooth solution (p, u) of (| 1 . 1 [) — (| 1 . 3p on R 2 x [0, T] for some 
time T > 0, with smooth initial data (po,uo) and smooth external force /, satisfying 

0<p<p (2.1) 

and 

A 1 +A 2 <2{C + C f ) e . (2.2) 
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In this paper, we choose 9 — i and assume that e < 1. 

Before proceeding, we remark that a careful application of the standard Rankine-Hugoniot condition 
to (jl.ip shows that discontinuities in p, P{p) and Vu across hypersurfaces can be expected to persist 
for all time, but that the functions F and w should be relatively smooth in positive time reflecting a 
cancellation of singularities (for example, see [BJ El El E] ) • We can rewrite the momentum equation in the 
form, 

pu } = djF + pd k w^ k + P p. (2.3) 

Thus L 2 estimates for pii, immediately imply L 2 hounds for VP and Vw. Stated differently, the decom- 
position (|2.3|) implies that 

AF = div(pu- pf). (2.4) 
These two relations (|2.3|) ~ (|2.4p will play the important role in this section. 



Proposition 2.1. There is a positive constant C = C(p) such that, if (p,u) is a smooth solution of 
£0])-£Op satisfying I2l]) - l2^) . then 

sup / \p\u\ 2 + G(p) dx + [ [ \Vu\ 2 dxdt < C(C + C { ). (2.5) 
0<t<TJ |_2 J Jo J 

Proof. Using the energy estimate, we can easily obtain (I2.5[) . and omit the details. □ 

The following lemma contains preliminary versions of L 2 bounds for Vu and pu. 

Lemma 2.1. If(p,u) is a smooth solution of i f 1. i)) ~i TO)) as in Provosition \2.1[ then there is a constant 
C = C(p) such that 



(2.6) 



sup a [ \Vu\ 2 dx + [ [ <rp\u\ 2 dxdt < C (C + C f + Ox) , 

0<t<T J Jo J 

where 0\ = f J a\Vu\ 3 dxdt, and 

sup a 2 I p\u\ 2 dx+ / / a 2 \Vu\ 2 dxdt 
o<t<T J Jo J 

< C , (C + C , / + A 1 (T)) + C^ J o- 2 {\u\ 4 + \Vu\ i )dxdt. (2.7) 
Proof. Multiplying (jl.lj) ? by aii, integrating it over R 2 x [0, t], we obtain 
J J ap\u\ 2 dxds 

= J J {—aii ■ VP + paAu ■ u + <rV((A + p)divu) ■ ii + apf ■ u) dxds 

4 

:= (2.8) 

i=l 

Using the integration by parts and Holder's inequality, we have 
Ji = cru ■ VPdxds 



-J J (-c(divu)t(P - P{p)) + a(u ■ Vu) ■ VP) dxds 
J adivu(P - P{p))dx 

II ( CTtdivu( - jP _ P ^ + <jPfdiYU + <J ( U ' Vm ) ' V - P ) dxds 
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= J adivu(P - P{p))dx - J J {a t divu(P - P{p)) 

+aP' (— p(divu) 2 — u ■ VpdivM + (u ■ Vu) • Vp) } dxds 
= J adivu(P - P(p))dx - J J {a t divu(P - P{p)) 

-aP'p(divu) 2 + aP(divu) 2 - aPd l w> dju'jdxds 

< C {u J |V<u||p - p\dx + J |Vujjp- p\dxds + J \Vu\ 2 dxds) , (2.9) 

J2 = J J apAu ■ iidxds 

= J \Wu\ 2 dx + ^ J \Vu\ 2 dxds- J ap l d i u j d i (u k d k u j )dxds 

< J \S7u\ 2 dx + ^ J J \Vu\ 2 dxds + CO l , (2.10) 



J3 = J J ctV((A + p)divu) ■ iidxds 



= -| J(X + p)\dWu\ 2 dx + j J p t \divu\ 2 + ^{\ + p)\d\vu\ 2 )dxds 

+ J J crV((A + /i)divu) • (u ■ Vu)dxds 
< J(\ + [i)\divu\ 2 dx + J ]^{\ + p)\divu\ 2 dxds + CO u (2.11) 



1 '* 



J a = I / apf ■ iidxds <- / ap\ii\ z dxds + C / a f ■'dxds, (2.12) 



where t G [0,T]. From ([27H]) . (|2~8| - ([2~T2| . we immediately obtain (|2~6|) . 

Takeing the operator 9 t + div(u-) in ![) ■?. multiplying by a 2 ii and integrating, we obtain 

a 2 f 

— j p\ii\ 2 dx 

+crV[6» i <9 t ((A + p)divu) + div(udj((X + p)divu))} + a 2 u 3 [{pf j ) t + div(upf r )]}dxds 
5 

i=l 

Using the integration by parts and Holder's inequality, we have 
K 2 = - J J a^ldjPt+dividjPu^dxds 
= J o J a2 \- d i ^ P'Pt + d k ii j dj Pu k ] dxds 

= J (j 2 [-P / pdivud j ii j + d k {d 3 u j u k )P - Pd 2 (d k u> u k )\dxds 
cQf J \Vu\ 2 dxds^j J a 2 \Vii\ 2 dxds^j 



< 
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< <7(C + C»*(/ / <7 2 \Vu\ 2 dxds) , (2.14) 



K 3 = J J pa 2 u j [Au J t + div(uAu j )]dxds 

= — J J a 2 \i\diV? diu\ + Av?u-Vv?]dxds 

a 2 p[\Vu\ 2 - d l u 3 u k d k d i u 3 - d i u j d i u k dkU j + Au 3 u ■ Vu 3 ]dxds 
a 2 n[\Vu\ 2 + diV? dku k div? - diU 3 diU k dkU J - diU 3 diU k dkU 3 \dxds 
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~\J Q J ^vlVufdxds + C J J a 2 \Vu\ 4 dxds, (2.r>i 



K4, = I I a 2 u J [djd t ((X + p)divu) + div(it<9j((A + p)divu))]dxds 



{a 2 dju 3 [d t ({\ + fijdivu) + div(u(A + /x)divu)] 
+a v? div(dju(\ + p)divu)}dxds 



a 2 djU 3 [dt({\ + p)divu) + u X'dkpdivu + (A + p)u k dkdivu]dxds + O4 

f { <7 2 djVp[{\ + p)——divu + X' p t divu + u k X 1 dkpdivu]dxds + O4 
Jo J ™ 

4 r D 

<j 2 dj{dtu J + u ■ Vm j )(A + p)-j^d\vudxds + O4 

cr2 ( A + M)l^ divu\ 2 dxds + Oi, (2.16) 

^5 = jf J <J 2 u 3 [(pf) t + div{upf 3 )}dxds 
= j J <T 2 u 3 [pfl + pu ■ Wf 3 }dxds 

< C f I [ap\u\ 2 + a 2 \u\ 4 + a 4 |V/| 4 + a 3 \f t \ 2 ]dxds, (2.17) 



where O4 denotes any term dominated by C f* J ct 2 |Vm| 2 (|Vu| + |-j=^divu|), and t € [0, T]. From (|2.13j) - 
(|2.17p . we immediately obtain (|2.7[) . □ 

The following lemmas will be applied to bound the higher order terms occurring on the right hand 
sides of rf2~6l)-d2~7l). 



Lemma 2.2. If (p,u) is a smooth solution of U.1)) - (T^\) as in Proposition ^. 11 then there is a constant 
C = C{p) such that, 

HlLP^C p (C + C7 / )||V«||^ a + Ci,(Cb + C , / )||V«||^, pe[2,oo), (2.18) 

||Vu||i, < C p {\\F\\ L v + \\w\\ LP + \\P - P(p)\\lp), V e (l,oo), (2.19) 
||VF|| iP + \\Vw\\ L v < C p (\\pu\\ LP + \\f\\ LP ), P e (l,oo), (2.20) 
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\\F\\ LP + \\w\\ LP < C p (\\pu\\J (||V W ||f 2 + ||P-P(p)||* 2 ) 

+\\Vu\\L* + \\fh* + \\P-P(p)h*), pe[2,oo). (2.21) 
Also, forO<t 1 <t 2 <T,p>2ands> 0, 

j* J a s \p-p\ p dxds<c(^J^ J a s \F\ p dxds + C + C^j . (2.22) 

Proof. Using the Galiardo-Nirenberg inequality, we have 

\\u\\ LP < C p \\u\\l 2 \\S/u\ff , pe [2,oo). (2.23) 

Since 

i 

p f \u\ 2 dx < f p\u\ 2 dx +( f |p - p\ 2 dx] 2 ( f \u\ 4 dx 



applying (|2.23|) and Proposition ^. 1[ we get 

IM| 2 L2 < C(C + Cf) + C(C + C»*||tt|| £a ||Vu|| ia 

and 

\\u\\h < C(Co + C f ) + C(C + C/)||Vu||| 2 . 

Combining it with (j2~2"3"]l . we get [[235]) . 

Since F = (A + 2p)div?i — P + P(p), we have 

,F + P-P (p) 
A + 2p 



Aul = d i( \ , o.. ) + d i( wJ ")- ( 2 - 24 ) 



From the standard elliptic theory we can get (|2.19p . 
We compute from (jl . that 

pAw^ k = d k (pui) - d 3 (pu k ) + djipf) - d k ( P P). 

Using the standard elliptic theory we can get 

IIVHUp < c p (\\pu\\ L p + \\f\\ LP ), pe(i,oo). 

From (12.31). we have 



djF = pv? - pd k w^ k - pf 

and 

||VF|| L p < C(\\pu\\ LP + \\Vw\\ L p + H/lli,), p € (l,oo). 

Thus, we can obtain (|2.20p . 

From (|2~20l) with p = 2 and (|2~23)) . we can immediately obtain (|2~2T|) . 
Multiplying the mass equation (|l.l|) i by pcr s \p — p| p sgn(p — p), we have 

(d t + div(u-))(<7 s |p - pT) + 3^((p - 1)P + p)| P - pMp - P(p)| 
= S a s -V t |p - p? - 3^((P - 1)P + P)sgn(p - p)|p - pr^. 
Integrating, we get 

*2 /-t 2 



fa s \p-p\ p dx +C'- 1 f 2 f a s \p- p\ p dxds 



<r(t 2 )Vti 



ti Jti 

-1 



(\J J a^\p-p\ p dxd S + J t 2 Ja°\F\ p dxd s y (2.25) 



Using the result of Proposition ^, li we can immediately get (|2.22|) . □ 
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Now, we apply the estimates of Lemma 12.21 to close the bounds in Lemma 12.11 

Proposition 2.2. If (p, u) is a smooth solution of il.l]) - FL3\) as in Proposition \2.1\ and e is small 
enough, then we have 



sup / (cr|Vw| 2 + cr 2 p\u\ 2 )dx + ( I ' (crp\u\ 2 + a 2 \Vu\ 2 ) dxdt <{C + C f ) 6 '. (2.26) 

<t<T J Jo J 



I (ct|Vm| 2 + cr 2 p\u\ 2 )dx + [ 

0<t<T 

Proof. From Proposition 12.11 and Lemmas 12. 1112. 2[ we have 

i-T 



LHS of (12211) < C(Co + C f ) + C J a J {o-\Vu\ 3 + a 2 \u\ 4 + cr 2 \Vu\ 4 )dxds. (2.27) 
From (|2.19p . we get 

J <T 2 \X7u\ 4 dxds < J a 2 [|F| 4 + M 4 + \P- Pip)]*] dxds. (2.28) 
From GTHfy and (|2~20|1 - (|2^3|1 . we obtain 

J J a 2 \F\ 4 dxds 

< C a 2 (J \F\ 2 dx^j (J \VF\ 2 dx\ ds 

< C sup [ o-{\X7u\ 2 + \p-p\ 2 )dx [ [ a (p\ii\ 2 + \f\ 2 ) dxds 

0<t<T J Jo J 

< C^+Co + Cj) 2 , (2.29) 

T r 

a 2 \w\ 4 dxds 



< C J a 2 (J \w\ 2 dx^ (^J \Vw\ 2 dx\ ds 

< C sup [ a\X7u\ 2 dx [ { a (p\u\ 2 + \f\ 2 ) dxds 

a<t<T J Jo J 



< 



C{Ai+C + C f Y, (2.30) 



jf J a 2 \p- p\ 4 dxds < c(j^ J a 2 \F\ A dxds + Co + Cf^ 

< C(A 1 +Co + C f ) 2 +C(Co+C f ), (2.31) 

<j 2 \u\ 4 dxds < C{C + C f )J cr 2 (||Vu||| 2 + ||Vw||i 2 ) 

< CA^Co + Cf) 2 + C(C + C f ) 2 . (2.32) 
From ![23g]H[2~32]) . we have 

T p 

a 2 {\u\ A + \Vu\ 4 )dxds < CA\ + C(C + C f ). (2.33) 



'0 



Similarly, we get 

J J cr\X7u\ 3 dxds < J J (cr 2 \X7u\ 4 + \X7u\ 2 ) dxds 
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< CA 2 + C(C Q + C f ). (2.34) 
Then, from |[2T2]) . (|2~27f and (|2~33l) - ([2~34|) . we obtain 

LHS of (l2~26l) < C(C +C f )+ CA\ 

< C(Co+C f )+C(C +Cf) 29 

< {Co + Cf) e , (2.35) 

when 

e^+e s <±. (2.36) 

□ 

Then, we consider the Holder continuity of u in the following lemma. 
Lemma 2.3. When t 6 (0,T] and a S (0, 1), we have 

< u(-,t) > a < C (||HI^(l|Vu||^ Q + (Co + Cf)^) + \\Vu\\ L * + (Co + C f )^) . (2.37) 

Proof. Let p — j^—. From (|2.19[) . (|2.2ip and Sobolev's embedding theorem, we have 

< u(-,t) > a 

< C||Vu|| £ , 

< C(\\F\\l P + \\w\\ Lp + \\P-P(p)\\ Lp ) 

< c (jHl5f (||Vu||f 2 + ||f - p(p)\\l) + \\Vu\\ L , + ll/iu, + \\p- ptf)\\v 
+c||p-p||| 2 || P -p||y 

< C MIHIJ 2 (||Vu||| a + (Co + C f )v) + \\Vu\\ L * + (Co + C f )*j . 

□ 

Proposition 2.3. If uq G H , (p,u) is a smooth solution of \1. 1\) - [T73\) as in Provosition \2.1{ then we 
have 

sup / \\7u\ 2 dx + I I p\ii\ 2 dxdt < C(l + M q ). (2.38) 

0<t<T J Jo J 

Proof. Using a similar argument as that in the proof of (|2.6|) . we have 

sup / \\7u\ 2 dx + [ [ p\u\ 2 dxdt < C(C + C f + M q ) + C [ [ \\7u\ 3 dxds. 

0<t<TJ Jo J Jo J 

Without loss of generality, assume that T > 1. From (|2.26p and (|2.34[) . we get 

sup / \X7u\ 2 dx+ [ [ p\u\ 2 dxdt<C(l + M q ) + C [ [ \X7u\ 3 dxds. 

<t<T J Jo J Jo J 



I \\7u\ 2 dx+ [ 

0<t<T 

From (pm?]) and (12~2"2"1) , we have 



sup [\Vu\ 2 dx+[ [ p\ii\ 2 dxdt <C(1 + M q )+C [ [( 

<t<T J Jo J Jo J 



0<t<T 

From (|2~2H|) - (|2T2T|) and ([2~2"3]) , we obtain 



cT t 

' (\F\ 3 + \w\ 3 )dxds. 



J (\F\ 3 + \w\ 3 )dx 
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< C^J \F\ 2 dx^j (^J \VF\ 2 dx^j 2 + C (^J \w\ 2 dx^j \Ww\ 2 dx^j ' 

< cj (|Vu| 2 + \p- p\ 2 )dx (J [p\u\ 2 + |/| 2 ) dx\ 2 . (2.39) 



Thus, from Proposition ^. 11 we have 



sup / \S/u\ 2 dx+ I I p\u\ 2 dxdt 
<t<T J Jo J 



0<t<T 

< C(l + M q )+C f 1 (\\Vu\\ 2 L2 +Co + C f )(\\^pu\\ L 2 + \\.f\\ L 2)d S 

Jo 

< C(l + M q ) + ^J^ J p\ii\ 2 dxdt + C WVuWizds. 



Using Gronwall's inequality, we obtain 
/ |Vw| 2 da;+ [ 

0<t<T 



sup / \X7u\ 2 dx + / / p\u\ 2 dxdt < C(l + M q )e% H v «ll^ ds < C(l + M q ). 

<t<T J Jo J 



□ 



Proposition 2.4. If uq £ H , (p,u) is a smooth solution of H 1. i)) -i TO)) as in Provosition \2.1[ then we 
have 

sup a / p\ii\ 2 dx+ / / a\Vu\ 2 dxdt < C(l + M q ). (2.40) 

0<t<T J JO J 

Proof. Using a similar argument as that in the proof of (12. 7|) . from (|2.38p . we have 

sup a f p\ii\ 2 dx + f f a\X7u\ 2 dxdt < C(l +M q )+C [ [a (|u| 4 + |Vw| 4 ) dxds. 

0<t<T J JO J Jo J 

Without loss of generality, assume that T > 1. From (|2.26p and (|2.33[) . we get 

sup a f p\ii\ 2 dx + [ [ a\\7u\ 2 dxdt < (7(1 + M q ) + C [ fa (|u| 4 + |Vu| 4 ) dxds. 

<t<T J Jo J Jo J 



0<t<T 

From (EU), (j2~T5)l - (|2~Tg|) . (|2~2"2")l and (|2~5g]) . we have 
cr / p\ii\ 2 dx + I 

0<t<T 

From (EH), (12~2TJ)) . (j2~23"l) and (j23gj) . we obtain 



sup cr / p|u| 2 dx+ / / a\Vii\ 2 dxdt < C{M q ) + C / / cr(|F| 4 + |w| 4 )da;ds. 

<t<T J Jo J Jo J 



J J a (|F| 4 + |w| 4 ) dxds 

< C J a (J \F\ 2 dx^ (j \S7F\ 2 dx^j ds + C J a (J \w\ 2 dx^j (^J \Vw\ 2 dx^j ds 

< cJ q ^(J (|V«| 2 + \P- P\ 2 )dx^j (J (p\u\ 2 + |/| 2 ) dx\ ds 

< C(l + M q ) + ||Vu||| 2 cr J p\ii\ 2 dxds. (2.41) 
Using Gronwall's inequality, we can finish the proof of this proposition. □ 
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Lemma 2.4. For any p S [2, oo), we have 

Ml* < C p ||Vpu||| 2 ||Vu||^ f + CpHVulU.. (2.42) 

Proof. Since 

p J \ii\ 2 dx < J p\ii\ 2 dx +f^J \p- p\ 2 dx^j 2 (J \u\ A dx^j 2 , 
applying (|2.23|) and Proposition 12.41 we get 

||u||| 2 <C||Vpu||| 2 +C||Vn||i 2 . 
From (|2.23[) . we can immediately obtain (|2.42p . □ 
Lemma 2.5. For any q € (0,2), we have 

i-T 



a—p\u\ 2+q dxds < C(l + Mq 2 ). (2.43) 
Proof. Let p = ■^z^- Using Holder's inequality, (|2.40[) and (|2.42[) . we have 



rT 



o- 2 2 9 p\u\ 2+q dxds 



I" 7 2 + „ 2(p-2-g) qp 

< C / a^||VP^|| i2 p - 2 M^ds 
Jo 

" ° Jo a ~^^^ 2 [\\Vpu\\U\Vu\\ L z p +||Vu|| l2 ) ds 

/ /-T \ 2 / rT \ 

c 



I <T||Vu||i a dt) (/ ||v^lli»*) I sup a\\^puf L2 ) 
Jo J \Jo J \te[o,T] J 

qp Ip — A — qp 

/ pT \ 2(p-2) / »T \ 2(p-2) / 

+C[ a\\Vu\\l 2 dt) / ||V^Hi»* sup <r|IVHl!= 



< C{l + Ml +i ). 



q2< YT^—^ V ^[0,p], (2.44) 



then we have 



sup a 

o<t<T J Jo 



' a 2+ i I p\u\ 2+q dx 



2 + q 



\^^-e 1+i o'p\u\ 2+q - (r 2+i \u\ q u 3 (djPt + div (djPu)) 



+pa 2+ i\u\ q u 3 (Au{ + div(uAu j )) 

+ct 2+ 'HW {d 3 d t {{\ + p)divu) + div(udj((X + ^)divu))) 



□ 



Proposition 2.5. If uq 6 if 1 , (p, u) is a smooth solution of il. i|) -i TO)) as m Proposition \2.4\ q <E (0,2) 



p\u\ 2+q dx + J cj 2+ i\u\ q \\Ju\ 2 dxdt<C{M q ). (2.45) 
Proof. Using a similar argument as that in the proof of (|2.13p . we have 
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+o 2+ i\u\ q u 1 {(pf)t + div{upf))}dxds 

5 

:= (2.46) 



i=l 

Using the integration by parts and Holder's inequality, we have 



H 2 = - J J a 2+ i {u^u 1 (djP t + div(djPu)) dxds 
= J a 2+ i\u\ q (d 3 u j P' Pt + d k ii J d 3 Pu k ) dxds 

+ Jo I a2+i P ' pt + d M q u 3 d 3 Pu k ) dxds 

= J q J ° 2+i {-\u\ q P'pdivudjU j + d k {\u\ q djU j u k )P - Pd j (\u\ q d k u j u k )) dxds 
+ Jo J (- d ^\ q u j P'pdivu + Pd k {d j \^ q ii j u k )- Pd 3 {d k \u\ q u j u k ))dxds 

= J J a 2+ i (-\u\ q P'd jU j pdivu + Pd k \u\ q d 3 u j u k + P\u\ q d j u j divu 

-Pd 3 \u\ q d k u j u k - P\u\ q d k u j d 3 u k - d 3 \u\ q u j P' pdivu 
+Pd k {u 1 u k )d ] \ii\ q - Pdj{u 3 u k )d k \u\ q ) dxds 

= J J ° 2+i ( - \u\ q P'd jU j pdivu + P\u\ q djU j divu - P\u\ q d k u 3 d 3 u k 

-d 3 \u\ q u j P' pdivu + PuPdivud 3 \u\ q - Pv? d 3 u k d k \u\ q )dxds 

< C j a 2+ i\\7u\\u\ q \\7u\dxds 7 (2.47) 

H 3 = J J p J a 2+ i\u\ q u 3 ' (Au\ ] +div(uAu : >) S J dxds 

= - J q J <? 2+i tJ- (di{\u\ q u j )diui + Au^u ■ V(|u|V)) dxds 
= - J <J 2+ h(\u\ q \Wu\ 2 - d t u^u k d k d t u 3 \u\ q - d l u 3 \u\ q d l u k d k u 3 
+d i \u\ q u>d i u> t + Au^u ■ V{\u\ q u j )} dxds 

<j 2+ ip{\u\ q \Vu\ 2 + diu j divudiu j \u\ q + di\u\ q u j divudiu j 



-d l u^\u\ q d l u k d k u ] + di\u\ q u j diu j - di\u\ q u j diU k d k u j 
-diU j diU k d k \u\ q u j - d l u j \u\ q d i u k d k u j )dxds 



< - J J a 2+ ip\u\ q \Vu\ 2 dxds- J* J l<j 2+ ip\u\ q - 2 \V\u\ 2 \ 2 dxds 

+C Jo J a2+ ^ Vu \ 2 ^ q \ Vil \ dxds ^ ( 2 ' 48 ) 

H 4 = J a 2+ ^\u\ q u : > (d 3 d t ((X + p)divu) + div{ud 3 ((X + p)divu)))dxds 
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{-a 2+ zdj{\u\ q u 3 ) {d t ((X + p)divu) + div(u(X + p)divu)) 
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ct 2+ 2 \u\ q u j div(d 3 u(X + p,)divu)}dxds 



-a 2+ id j (\u\ q u j )(d t ((X + p)divu) 



+u X'dkpdivu + (A + p)u k dkdivuj dxds + O5 

J J -a 2+ ^dj(\u\ q u j ) ((X + p)^divu + X'p t divu + u k X'd k pdivu^j dxds + 5 
- J J <J 2+ ^\u\ q dj(d t u ] + u ■ Vu J )(A + fi)^divudxds + 5 



D 



a 2+ ^dj\u\ q vP{X + n) — divudxds 
.'0 J Dt 

2 



< 



2+f 



(A + p)\u\ q 



+ J / qa 2+ HX + p)\u\ q 



— - divu 
Dt 

— - divu 
Dt 



dxds + O5 
\X7ii\dxds, 



(2.49) 



cr 2+ i\u\ q u J [(pf r ) t + div(upf )}dxds 



a 2+ i \u\ q u j [pf t + pu ■ Vf j ]dxds 



C 



|2+9 



\ft\ 



2+9 1 U,|4+2g 



|V/| 4+2 «) dxds, 



(2.50) 



where O5 denotes any term dominated by C f Q J <t 2+ 2 Vit| 2 (| Vu| + | jjj divu \)dx ds, and i € [0,T] 
From gag-QHS), (l2~2TD - (l2"22l . (12351) . (MB, (gl2| - (|13g) and (lAl - (l23(JD . we have 



sup a 2+ 5 / p|u| a+ «da; + / / <j 2+ %\ii\ q \X7 u \ 2 dxdt 

0<t<T J JO J 

< C(M q ) +C J J a 2+ i (\ii\ q \X7u\ 4 + \u\ q \Vu\ 2 + \u\ A+2q ) dxds 



< C{M q )+C 



I 2-9 

Il 4 +< 



+C J (\\Vuf L t 2 ' 1 + \\Vut+ 2,! ) 





V 2 


2 


- q ds 1 




9(2-9) 




\ 4(2 + <j) 


(is 




2g 


+ ||Vu| 









r- 




^ (711^1^4+9^1 









(2-9)(4+9) 
4(2 + 9) 



(j^||Vm||[V <!)(4+5) rfs 



/•T / _j_ 2(2 + 9) \ 

< C(M q )+Cj^ a^W^uW^PWVuW^ + \\Vu\\ 2 L2 jdt 

f T _J_ 

+C I a 2 -" (\\F\\ L 4+ q + \\w\\ L i+ q + \\p- p\\ L i+ q ) 2 -> ds 
Jo 

< C(M q )+C I a^{\\F\\ L , +q + \\w\\ Li+q )^ds + C I a^\\p - 



ds 







< C(M q )+Cj a— (||F||£ 9 ||VF||£« 
1 ., 

~7=77 II Fl, 

li 4 +9 



L 2 



I Viol 



2 + 9 \ 2-9 



£2 



ds 



/ T (T^||F|| r +l& 
Jo 
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f T 4 16 8(2 + .;) 

< C{M q )+C / <7*=«(||Vu||l» + pH^a) l«+«Ha-« (Hpulli, + ||/|| £a )l«+«H»-«ti8 

Jo 

+c / ^(||Vu|U 2 + ||p-p|U 2 ) 2 (||pti|U 2 + ||/|| L2 ) 2 +«rf s 

Jo 

< C(M g ). 

□ 

Proposition 2.6. If f E L^°L^ q , (p,u) is a smooth solution of 1 1. Jj )-i TOp as in Proposition \2.5\ then 
we have 

\\F\\ L oe + \\ W \\ L oo < C(\\Vu\\ L 2 +\\ P - p\\ L 2)^(\\p U \\ L 2 + , + \\f\\ L 2 + q)^ (2.51) 

and 

(||F|| L » + \\w\\ L -)da < C(M q )(C + C/)*f (1 + T). (2.52) 
Proof. From f2~5]l . (|2T2U|) . ([2~2^1) . (j2~4l))) and the Galiardo-Nirenberg inequality, we have 

< C||F||^||VF||^ 

g 2+q 

< C(|| Vu|| L 2 + \\p - p\\ L ^ 2 + 2 " (\\pu\\ L 2+ q + 2 + 2 " ■ 



and 



/ \\F\\ L ~ds 
Jq 



i-T 

11 6 q 4 + q 2 + q 

< C(M q ) (CT _ 2(C + C , / )2)2 + 2,( cr 4+2, ) 2 + 2, rfg 

Jq 

q6 2 + q qe 

< C(M q )(C + Cf)— <j-—*ds< C{M q )(Co + C f ) — (l + T). 

Jo 

Similarly, we can obtain the same estimates for w. □ 

Then, we derive a priori pointwise bounds for the density p. 

Proposition 2.7. Given numbers < p 2 < p 1 < p < p\ < p2 < p, there is an e > such that, if {p, u) 
is a smooth solution of U. ij )-) TO)) with Cq + Cf < e and < po < pi, then 

0<p<p 2 , (x,t) e M 2 x [0,T\. (2.53) 

Similarly, if po > p^ for all x, then p > p 2 for all x and t. Furthermore, the estimates in Propositions 
dUtH hold. 

Proof. At first, we prove that if (|2.1[) and (|2.2[) hold, then estimate (|2 . 53[) holds. 
We fix a curve x(t) satisfying x = u(x(t),t) and x(0) — x. From (ll.21|) . we have 

j t A(p(x(t),t)) + P(p(x(t), t)) - P(p) = -F(x(t),t), (2.54) 

where A satisfies that A(l) = and A'(p) = 2m + a(p) . 

First for small time, we estimate the pointwise bounds of the density as follows. From (|2.ip and 
lf2"32]> . we have, for all t e [0, 1], 

\A(p(x(t),t))-A(po(x))\ < C(M q )(Co + C f )^ +Ct. 

When 

qff | 

2C(M q )e—* < A(pi + -(P2- Pi)) - A(pi), (2.55) 
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we get 

A(p(a(t), t)) < A(pi + ~(p 2 - p x )), t e [0, r], 

and 

p{x,t) <pi + -(p 2 -pi), (x,t) £t 2 x [0,r], (2.56) 
where r = min{l, 2^(A(pi + |(p2 — pi)) — A(pi)). Since po > 0, then we have 

A(p(x(t),t)) > A( Po (x)) - C(M q )(C + C/)*& - Ct > -oo, i e [0,r], 

and 

p > 0, (x,t) £l 2 x [0,t]. 

Similarly, if po > p l and 



we get 



2C(Af 9 ) £ ^ < A(pJ - A( El - i(p x -p 2 )), (2.57) 
-°^£i-^(£i-£2)' (^,*)6K 2 x[0-n]. (2.58) 



3 

where n = min{r, ^(AQ^) - A(p x - i(^ - p 2 )))}. 

Then, for large time t > r, we estimate the pointwise bounds of density as follows. From (|2.5j) . (|2.26[) . 
(j2~45ll . ([23T|) and $TEfy . we have 

d -^^3 + P(p( x (t),t)) - P(p-) = 5 (t), (2.59) 

where 

\O 5 (t)\<C(T,M q )(C + Cf)^, t>r. 
Now, we apply a standard maximum principle argument to estimate the upper bounds of density. Let 

to = max{f e (t,T}\A( P (x(s),s)) < A(p 2 ), for all s € [0,t]}. 

If to < T, we have 

A(p(a;(to),to)) = A(p 2 ), 
dA(p{x(t),t)) 



>0, 

t=t 



and 

From (|2.59p . we have 
On the other hand, when 
we have 



dt 

p(x(t ),t ) = p~2- 
O 5 (t ) > P(f) 2 ) - P(f>). 



C(t, M q )e^k < P(p 2 ) - P(p), (2.60) 



O 5 (t ) < P(p 2 ) - P(p). 
It is a contradiction. Thus, we have to = T and 

P < P2, £l 2 x [0,T]. (2.61) 

Similarly, we can obtain the lower bound of the density. 

Using the classical continuation method, (|2.26p and (|2.53p , we can finish the proof of this proposition. 

□ 

Then, we can prove the global existence of smooth solutions to (|l.ip - (|1.3[) . 
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Proposition 2.8. Assume that p Q - p S W 1 ^ D C 1+a , u E H 2 n C 2+a , p > 2, a £ (0,1), poO) > P 
with some p > for all x e R 2 , P, A G C°°([0, p]) and / S C°°([0, oo); C°°). Under the assumptions 
of Theorem XTJ\ i/ien tfiere exists a solution (p,u) e C 1+Q! ' 1+Q / 2 x C 2 +" 4+q / 2 (M 2 x [0,T]) satisfying 
\1. 1\) - H77S\) and for which the bound estimates of Provositions \2. hold, for all T > 0. 

Proof. Using similar arguments as that in [16j and in the proof of Proposition 3.2 in [8], one can obtain 
this proposition. □ 



3 Proof of Theorem 11.11 

Let js(x) be a standard mollifying kernel of width 5. Define the approximate initial data (pq, Uq) by 

Po = 3s * Po + 8, u 5 o = 36 * Mo- 
Assuming that similar smooth approximations have been constructed for functions P, f and A, we may 
then apply Proposition 12.81 to obtain a global smooth solution (p s ,u s ) of (|l.l|) - (jl.3|) with the initial data 
(p^, Uq), satisfying the bound estimates of Propositions 12. lj | 2. 61 hold with constants independent of S. 
First, we obtain the strong limit of {u s }. From (|2 . 26[) and (|2.37[) . we have 



< «*(•,*) > a < C(r), f>r>0, ae (0,1). 



(3.1) 



From (13.11). we have 



u s {x,t) - 



Br(x)\ 



Br(x) 



u s (y,t)dy 



< C(r)R a , t>r>0. 



Taking R = 1, from f2~T5|) and (|2~2^|) . we have 



\\u ||l«=(R 2 x[t,oo)) < C(t). 

Then, we need only to derive a modulus of Holder continuity in time. For all t2 > t\ > r, from 
(j2Tl2l) and ((221), we have 



(3.2) 



< 



\u d (x,t 2 )-u d (x,h)\ 

1 



\B R (x)\ J tl J Br {x) 



\u 5 t (y,s)\dyds + C(T)R a 



< CR-^h-t^ 



u s t \ 2 dyds) +C(t)R° 



< CR^\t 2 - ti.fi (J 2 J \u 5 \ 2 + \u s ■ Wu 5 \ 2 dyds 

< Cir^R-^h-hfi +R a ). 



C{T)R a 



Choosing R — \t% — ti \ 2 + 2 ° , we have 

<«' >^oo)<C{r), r>0. 
From the Ascoli-Arzela theorem, we have (extract a subsequence) 

u s -> it, uniformly on compact sets in R x (0,oo). 



(3.3) 



(3.4) 



Second, we obtain the strong limits of {F 5 } and {w 5 }. From ([2~20]) - ((2~2T|l . (|2~26| and ([2~45ll . using 
similar arguments as that in the proof of (|3.ip ~ ()3.2p . we have 



< F 5 (-,t) > a +||F 5 || ioo(K2x[T , T]) + < w 5 (-,t) > a +|h/|| L =c (R 2 x[T , T]) < C(t,T), 



(3.5) 



20 



where < r < t < T and a' G (0, 5+^]- The simple computation implies that 



F* = p d (2p + \(p*))\F 



« d 



rfs \2/x + A(s) 



+ d_ / P(s)-P(p) 
s=pS ds \ 2 M + A(s) 



divu 5 



-u* • VF S + (2/i + A(p 5 ))divw 5 - (2/i + A(p 5 ))a i M^ i wf (3.6) 

and 

(w 5 ) k t > j = -u s ■ V{w 5 ) k ' j + djiii - d k u] - djuldiut + d k utd lU 5 r (3.7) 
Then, from (j2~20| . (|2~26j) . (|2~33| . (f3T2|) and ([33]l . we have 

ll-Fj?llz 2 ()8 2 x [t,t]) + II IU 2 (r 2 x[t,t]) < C(t,T), T > t > 0. 
Using a similar argument as that in the proof of (|3.3j) . we obtain 



< > K ^] + < ™* > R ^]< C(r,T), T > r > 0. (3.8) 

and (extract a subsequence) 

— > F, w s — » u>, uniformly on compact sets in R 2 x (0, 00). (3.9) 

Third, we obtain the strong limit of {p 5 }. From (|2.53p . we get (extract a subsequence) 

p A ' A p, weak-* in L°°(R 2 ). 

Let &(s) be an arbitrary continuous function on [0, p]. Then, we have that (extract a subsequence) §(p s ) 
converges weak-* in L°°(R 2 ). Denote the weak-* limit by $: 

#, weak-* in L°°(R 2 ). 

From the definition of F, we have 

divu = PF + Po, (3.10) 

where 



"^ = XT? V P °to = "(pX^OO - TO)- 
2/i + A(p) 



From Ijl.ip . we have 



dtphip + div(p In pu) + Fpv + pPo = 

and 

d t (p In p) + div(p In pu) + Fpv + pP = 0. 



Letting ^ = p In p — p In p > 0, we obtain 

9** + +div(#n) + F{pU - pv) + Fp(i/ - p) + pfb - pH = 0. (3.11) 
with the initial condition ^^=0 = almost everywhere in R 2 . Let 4>(s) = sins. Since 

^'(«) = -> *€[0,p], 
s p 

we get 

0(p 5 ) - 0(p) = 0'(p)(p 5 - p) + 1 -4>"{ P + t(p s pW p) 2 ' £ g [0, i], 



and 



lim||p*-p|| 2 2 <C||*|| L1 . (3.12) 
a— >u 
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Similarly, every function / £ C 2 ([0,p\) satisfies 

.</(/ f(p))dx <cj \g\*dx, 
where g is any function such that the integrations exist. Then, when v £ C 2 ([0,/5]), we have 

< C ( \F\^dx 



and 



When P e C 2 ([0,p]), we have 



F{pv — pv)dx 
Fp(y — v)dx 



< C / \F\Vdx. 



(pPo - pPa)dx 



< 



(pP - pPo)dx 



p(P - P )dx 



< C / ^dx. 



(3.13) 

(3.14) 
(3.15) 

(3.16) 
(3.17) 



When P is monotone function on [0,p], using the Lemma 5 in |16] , we have 
From ((3~TT|) - (l3~T7l) . we obtain 

J ®dx < J {1 + \F\)^dxds. 
Using (|2.52p and Gronwall's inequality, we get 

* = 0, (t,x) S [0,T] x M 2 , 

and (extract a subsequence) 

p s - p — > p - p, strongly in L k (R 2 x [0, oo)), 

for all k e [2,oo). 

Thus, It is easy to show that the limit function (p, u) are indeed a weak solution of the system 
HH)-([I3]). This finishes the proof of Theorem O □ 



4 Lagrangean Structure 

Proof of Theorem [TT2l 

(1). Here, we consider the case to = 0. The proof of the case t > is similar, and omit the details. 
We denote that X(t, x ) = X(t; x , 0). 

To prove the existence of the integral curve X(-, xo), we first assume that X s (-, xo) is the corresponding 
integral curve of u s with initial point xo el 2 , 

X s (t,x ) =x + [ u s (X 5 (s, x ), s)ds, se [0,oo). (4.1) 
Jo 

From (|2.18p . (|2.3T[) — (|2.38[) and (|2.40|) . using a similar argument as that in the proof of (|3.2p . we have 

< u s (-,t) > a < C + Ca-%, a e (0,1), (4.2) 

||<AV)IU» <C + Ca-f, (4.3) 

and 

/ \\u s (t, OHi-dt < C(T -Tx + T 1 -^ - Tl ^), < Ti < T < oo. (4.4) 
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This implies that X 5 (-,xq) is Holder continuous on [0, oo), uniformly in S. Therefore, there is a subse- 
quence X Sj (-,xo) and a Holder continuous map X(-,xo) such that X Sj (-,xo) — > X(-,xo) uniformly on 
compacts sets in [0,oo). From this uniform convergence and (|3.4p . we have that X(-,xo) satisfies ()1 .22|) . 

Next we prove the uniqueness and continuous dependence for integral curves of u. Thus, let Xi(-, yi) 
and X 2 (-, y 2 ) be any two integral curves of u with respective initial points yi,y 2 € K 2 and define 

\u(X 2 {t,y 2 ),t)-u(X 1 {t,y 1 ),t)\ 

git) = 



and 



for t E [0, 00), where 



m(\X 2 (t,y 2 )-X 1 (t,y 1 )\) ' 
u s (X 2 (t,y 2 ) 1 t)-u s (X 1 (t,y 1 ),t)\ 



m{\X 2 {t,y 2 ) - X^y^) 

1 x(l — Inx), < x < 1, 
mix) = < ' ^ 

' x, 1 < a; < 00. 



From (fL22|) . we get 

\X 2 (t,y 2 )-X 1 {t,y 1 )\<\y 2 -y 1 \+ f g(s)m(\X 2 (s,y 2 ) - X 1 (s,y 1 )\)ds. (4.5) 

Jo 

Denote by < v >ll the log-Lipschtz seminorm of a given vector field v on R 2 : 

\v(x) — v(y)\ 
<v> LL = sup J-i> 

x=£y 

From (|2.24p and Proposition 2.3.7 in [3], we have 

g S (t) << u 5 (-,t) > LL < CH^IIfl^^ < C(\\u s \\ L 2 + \\F S \\ L ~ + \\p 5 - p\\ L - + \\w 5 \\ L ~). 
From |Ll3|) . (j2~5|l . (|2~T8l) and ([2~52]) . we obtain 

/ £M<C(1+T). 
Jo 

Using (|3.4p and Fatou's lemma, we have 

l(s)ds < liminf / g s (s)ds < (7(1 + T). 



(5^0 

Using Gronwall inequality in (|4.5p . we get 

\X 2 {t,y 2 ) - X^yx)] < cxp(l - e~ K gds )\y 2 - y x \^- & ^ 

< eMi-e- c{1+T) )\y2-yir pi ~ c(1+T)) - (4.6) 

Thus, if yi = y 2 , then X\ = X 2 . There is therefore exactly one integral curve originating from a given 
point at time t = 0. From this uniqueness, we obtain that the convergence X s (t, yi) — > X(t, yi) uniformly 
on compact sets in [0,oo) for entire sequence 5 — > 0, independently of y±. 

(2). We prove the injection of X at first. Suppose that X(t, j/i) = X(t,y 2 ) for some t > and 
Vi, V2 S K 2 . Then for any ,s G [0, t), 

\X(s, yi ) - X(s, y 2 )\ < J < u(r, •) >ll m(|X(r, yi ) - X(r, y 2 )\)dr. 
Using a similar argument as that in the proof of (|4.6| . we have that X(s,yi) = X(s,y 2 ) for all s£ [0, t]. 
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Next we show that X(t, -)\q is an open mapping. Let A be an open subset of Q, y\ €. A and 
B ri {yi) C A, Q < s < t, zi — X(t,yi). From (|1.22p . using a similar argument as that in the proof of 
(|4.6|) . we have 

\X(s, yi) - X(s; z,t)\< exp(l - e^ 1 ^)^ - z |^p(-C(1+*))_ 
Thus, there is a sufficient small constant r 2 such that, if z £ B r2 (z\), then 

\yi-X(0;z,t)\< ri . (4.7) 

Thus X(0; z, t) £ -B ri (yi) C A if z e i? r2 (zi). From the uniqueness of the integral curves of u, we obtain 
that z — X(t, X(Q; z, t)) £ -)A Therefore, B r2 (z\) C X(t, -)A, and X(i, -)|o is an open mapping. 

(3) . Using a similar argument as that in the proof of (|4.6|) . we have 

\X(t 2 ,y 2 )-X(t 2 , yi )\ < exp(l - e~^° ds )\X{t llV2 ) - X^y^- ^ ^ 
< expll-e-^ 1 ^))^^,^)-^^,^)!^-^ 1 ^, 

which proves part (3). 

(4) . Part (4) of Theorem If .21 is an immediate consequence of part (3). □ 

To prove Theorem 1 1.3[ we need the following lemma. 

Lemma 4.1. Given x G X(t, -)V , say x = X(t,y) with y £ V, there is a sequence 5j — > 7 which may 
depend on x, such that X Sj (t, ■)~ 1 (x) tends to y as Sj — * 0. 

Proof. Using the argument as that in the proof of Part (1) of Theorem 11.21 we have that integral curves 
X s (s; x, t) of the approximate velocity field u s , defined by 

X s (s;x,t) = x- J u s (X s (r;x,t),T)dT, 

are Holder continuous in s e [0,t], uniformly with respect to 6. Therefore, there is a sequence Sj — > 
and a map X £ C([0,£];R 2 ) such that X j (-;x,t) converges uniformly to X(-), which satisfies 

X(s)=x- ( u{X{r),T)dT, 0<s<t. 

J s 

From the uniqueness of integral curves proved in Part (1) of Theorem II .21 we have that X(t) = X(t,y). 
Taking s = 0, from (|3.4|) and (|4.3j) . we get that y Sj — X Sj (t, ■)~ 1 (x) converges to y as Sj tends to zero. □ 

Proof of Theorem [TT3l 

Applying a standard maximum principle to the mass equation, using a similar argument as that in 
the proof of Proposition 12. 7\ we have 

p 5 (X 5 (t )2 ;),i)>p-, 

for any y S V, and all 5 sufficiently small. Let x — X(t, y) S V* and y s — X s (t 7 ■)~ 1 (x). From Lemma 
14. 1[ we have that there is a sequence Sj tending to zero such that y Sj tends to y S V. Then, for all 
sufficient small Sj, we get 

y s ' G V and ft (x, t) = p 5 i (X^ (t, y 5 >),t) > p~. 
From the convergence p s — > p (which holds for a. a. x), we obtain that 

p(x,t)=p(X(t,y),t)>p-, 

for all y £ V. From Part (2) of Theorem 11.21 we can finish the proof of this theorem. □ 

Proof of Theorem 11.41 

For any y £ U, there is a sufficient small So such that 

dist(y,dU) > 2S a . 
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Let Ug = {x G U\ dist(x, dU) > So}. Then, we have 

P S o\u So =S, V5<S . 

From (pH>2")) and (12"34"|) , we have 

A(p s (X s (t,z),t)) < A(p s a (z))- f P(p s (X s (s,z), S ))ds + tP(p)+ f \\F s (;s)\\ L ~ds 

Jo Jo 
< A(S)+C{T), 

for all te[0,T],S<5 and all z G U 5(> . Since A(C5) - A(S) > 2/x(ln(C5) - In J) = 2^1nC*, then we have 
btf(X% z), t)) < A(5) + A(C(2» - A(S) = A(C(T)5), 

and 

pP{X%z),t)<C(T)8<C(T)8o, 

for all t G [0,T], 5 < S and all z G U S(r Let a; = X(t,y) G U l and y s = X s (t, -^(x). From LemmalHJ 
we have that there is a sequence Sj tending to zero such that y Sj tends to y G Ug . Then, for all sufficient 
small Sj, we get 

y s > G t/ 5o and = p*{X s >(t,if>),t) < C(T)5 , 

for all t G [0, T]. From the convergence p s — > p (which holds for a. a. x), we obtain that 

p(M)=/oTO,y),t)<C'(T) ( 5o, 

for all t G [0, T]. Letting 6 -> 0, we get that j(?(X(t, y), t) = for all i G [0, T]. From Part (2) of Theorem 
O we obtain p(-,t)|[/' = 0, t G [0,T]. 

From flUggD , (|2T2%|) and (|2~2^|) , we have 



; (|p 5 -/5| 4 + |F 5 | 4 )^<C. 
From the convergence of {p s } and {F s }, we get 

J™ j(\p-p\ 4 + \F\ A )dxdt<C. (4. 

From (|2.25p . we have 

/> /> /-iV+l /> 

y |/-p| 4 (a:,f)da:< y |/ - p| 4 (a;, s)dx + C J J \F s \ 4 dxdt, 
where s,t € [N, N + 1], N > 1. Integrating it with s in [N, N + 1], we obtain 

>.JV+1 



//■JV+l /■ 
|// - p| 4 (x, t)d!B < C / / (|/ - p| 4 + |F 5 | 4 )dX(it 

From the convergence of {p s } and {-F 5 }, we get 

sup / \p - p\ 4 {x, t)dx <C {\p-p\ i + \F\ A )dxdt. 

Letting N — > oo, using (|48| . we can easily obtain (|1.24|) - (|1.25|l . □ 
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5 Propagation of Singularities 

Before proving Theorem 1 1.5i using the similar method in [6], we give the following three lemmas. 

Lemma 5.1. Given xq € K 2 and R > 0, there are positive constants 5q and ro, and a subsequence 
S = Sj -> 0, such that X 5 (t, ■)- 1 B ro (x t ) C B R (x ) for all S < S . 

Proof. By Lemma [4.11 there is a sequence 5^0 such that j/q := X s (t, ■)~ 1 (x t ) tends to Xo as S — > 0. 
Therefore, it suffices to show that there is a positive constant ro such that \X 5 (t, ■)~ 1 (x) — t/q | < R for 
sufficiently small 5 and for all x € B ro (x t ). 

Letting y s = X 5 (t, from ijiTTj) . we have 

X s (s;x,t)-X s (s;x t ,t) = x- [ u s (X s (r; x,t),r)dr - x* + / u 5 (X 5 (r; a; ,£), r)dr, 



for any s G [0, i). Using a similar argument as that in the proof (|4.7[) . we have 

|2/ 5 - <exp(l-e^ 1+t ))|x 
of sufficiently small radius rg. Then, we finish the proof. 



c ( 1+t >)|x* - x \°M-c(i+t)) <R 



□ 



Lemma 5.2. Let 5 — > &e i/ie sequence of Lemma \5.1\ then for all ro > 0, •) 1 — » •) 1 

uniformly on B ro (x ). 

Proof. Let y s := O -1 ^) and y := X(t, O" 1 ^) for x G S ro (x ). For < s < we have 

|X 5 ( S ,y 5 )-X( Sl y)| 

u s (X s (T,y s ),T)dT-x + / u(X(r,y),r)dr 



/ | M 5 (X A "(r,y 5 ),T)- U 5 (X(r,y),T)Mr + ^ \u(X(r, y), r) - M 5 (X(r, y), r)\dr 
< [ g S (r)m(\X s (r,y s ) - X(r,y)\)dr + [ \u(X(t, V ),t) - u s (X(r,y),r)\dT. 



Using a similar argument as that in the proof (|4.6[) . we obtain 
\y s -y\ = \X s (0,y s )-X(0,y)\ 

< exp(l- e - c ( 1+T ))^ \u(X(r,y),T)-u 5 (X(r,y),r)\dA 



exp(-C(l+T)) 



From (|4.4p and the uniform convergence of u 5 to u on compact sets in K 2 x (0,oo), we can get that 
y S (x) — > 2/(2;) uniformly on _B ro (x ). □ 

Lemma 5.3. Given ro sufficiently small and given t > 0, t/iere is a nondecreasing function 77 : [0, 00) — * R 
satisfying lim r ^o = such that, for 5 as in Lemma \5.2\ sufficient small, 



\X s (t 2 ,y 2 ) - X s (t 2 , yi )\ < V (\X s (t 1 ,y 2 ) - X 5 (t X , 
where h,t 2 <E [0,t] and X s (t,yi), X s (t,y 2 ) E B ro (x t ). 
Proof. We have that, for < ti, t 2 < t, 

\X s (t 2 ,y 2 )-X s (t 2 , yi )\ 



< \X s (t 1 ,y 2 )~X\t 1 ,y 1 ) 



\u s (X s (s,y d 1 ),s)-u d (X d (s,y" 2 ),s)\ds 



81 v 8i 



< \X s (t 1 ,y 2 )-X s (t 1 ,y 1 )\+ /( S )m(|X 5 ( S ,^)- X d (s,y 5 2 )\)ds. 

Jo 
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Using a similar argument as that in the proof (|4.6p . we obtain 

\X s (t 2 , y 2 ) - X s (t 2 , Vl )\ < exp(l - e- c ^)\X s (h, y 2 ) - X s (h, J/1 )| oxp( - c(1+t)) . 
This finishes the proof of this lemma. □ 

Proof of Theorem 11.51 

First, we show that 

osc(p; xl, E l ) < osc(p; x£, (E n B) v ) = 0, (5.1) 

for some open ball B centered at xq, where (EnB) 1 :— X(t, -)(EC)B). In deed, since the map X(t, : 
V* — > V is continuous by Theorem 11.31 (2), there is a positive constant r B such that B ri (x\ ) ) C B* for 
< ri < r B . This implies that J5* n B ri (x ) c(Efl £)*, so that, for any R > and ri < min{i?, r B }, 

(esssup - essinf)p(-,t)| B t nBri ( a .i) < (esssup - essmf)p(-,i)|( BnB )t nBfl(x *). 



Letting first rj — > and then R — > 0, we can obtain (|5.1|) . 
Case 1. If lim^_,o essinfpo| B n.B R (xo) = then for any So > 0, there exists tq > such that po|,EnB ro ( 
Sq. Using a similar argument as that in the proof of Theorem 1 1.4i we have 

p\{EnB T(l {x )Y < C(T)S , 

and 

osc(p{;t);x t ,(EnB ro (x )y)<C(T)6 , te (0,T]. 

From (|5.ip . we have 



jo 



< 



osc(p(-,t);4,£: t ) < C(T)S , < t <T. 
Letting So — ► 0, we get 

osc(p(-,t);4, J E*) = 0. 

Case 2. If lim^^o e ssinfpo|Bns n (a; ) > 0, then there exist ro > and p > such that po|_Ens rf) (x ) — P- 
Let -B be an open ball centered at xq such that B C B ro (xo). By Theorem[L2](2) and Lemma [BTTl we 
can choose a positive constant r% G (0,ro) such that, if Xj = X(t,ijj) = X s (t,yj) e B ri (x ) C\ (E (1 -B)*, 
j = 1, 2, then j/,, j/| S B n £7. From the mass equation (|1.1|) for (p 5 , u 5 ), we obtain 

^A(/(X 5 (i; -,0), t))\H = -P(p S (X s ( t] •, 0), t))|£ - F 5 (X 5 (t; 0), t)|£ 



= -a(t)A(p s (X s (t;;0),t))\y y l - F s (X s (t; -,0),^ 

_ P(p S (X s (t;-,0),t))\y y 2 



where a(t) = , S , YS ,' \^\ xu vl ■ From (|1 . 13[) and Theorem 11.31 we have \a(t)\ < C. Using Gronwall's 
inequality, we get 

\A(p s (x 2 ,t))-A(p s ( Xl ,t))\ (5.2) 
< C(T)\A(p s (y 5 2 )) - A(p s M))\ + C(T) f \F S (X S ( S , y s 2 ), s) - F s (X s (s, £\a)\d» 



< C(T)|A(p*(y 2 5 ))-A(^(^))|+C(T) sup \X s (s, y s 2 ) - X s {s, y{ )|*f« [ g F .da, 

0<s<t Jo 

where gF s i s ) =< F s (-, s) >^t q ■ Using a similar argument as that in Proposition 12.61 we have 

/ g F sds < C{T). (5.3) 
Jo 

From Lemma 15.31 we obtain 

sup |X 5 ( S ,y!)-X*( S ,yf)|^ < ( v (\x 2 - Xl \))^ < fa(2ri ))»*«. (5.4) 

0<s<t 
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Since yj G E n B, there is a r 2 > such that B r2 (y 3 ) d E H B. By Lemma [531 y| G Bra (%) for 5 
sufficiently small. Thus, if <S < ^, then |y - < |y - y|| + |y| - yj\ < S + ^ < r2 for all y G B s (yj); 
that is B<y(yf) C Br 2 (yj) C E B. Also, by lemmas for y G £<5(yf), 

|y - a?o| < |y - i#| + \y- -y s \ + \yj -x \<5 + \ y s 3 - Vj \ + »,(n) < 2» ? (n) ) 

that is, Bs(yj) C B2^( ri )(xo) f° r ^ sufficiently small. Thus, Bs(yj) C (£ fl B) fl B 2ri { ri ) [xq)- Then, since 

p s o(y S j) = $B S (yD Mvl ~ y)p 5 o(y) d y + 5 , we have 

essinfpo|(Bnfl)nfl 2 , (ri) (x ) < PoiVj) ~ $ < esssupp 1 (£nB)nB 2 „ (ri )(*<,)■ (5-5) 

From (|5.2 p -([5T5] l . we obtain 

\A(p s (x 2 ,t))-A(p s (x 1 ,t))\ < C(T)(ess S up-e SS mi)p Q \ {EnB)nB2Mri)(xo) + C(T)r,( ri )^ , 

for x\,X2 G [Ef] BY n B ri (xq), sufficiently small r\ > and sufficiently small <5 > 0. Taking the limit as 
6 — > 0, we get 

|A(p(x 2 ,t)) - A(p(si,t))| < C(T) (esssup-essinf)^ | (snB)nB^ (ri) (x ) +C(T)»j(ri)5fc, 

for xi, x 2 G (i?n B) f n B ri (xq) and sufficiently small r\ > 0. Then, we have 

(esssup-essinf)p| (BnB)tnBri(a .t) < C(T) (esssup - essmf)po|( EnB)nB2i)(ri)(xo) + C(T)r?(ri)s+5, 

for sufficiently small r\ > 0. Letting ri — > 0, using the condition osc(po; xoj E Hi?) =0, we obtain 
osc(p(-, t); x , (E n BY) — 0. From (|5.ip . we complete the proof of Theorem 1 1.51 □ 

Proof of Theorem 11.61 

It follows immediately from Theorem 11.51 that . under the conditions of Theorem 11.61 (a). p(-,t) has a 
one-sided limit with respect to M 1 :— X(t, at the point X(t, xq). Then, since F = (A + ^i)divw — 
P(p) + P{p) is Holder continuous for t > by (|1.18p . divu(-,£) has a one-sided limit with respect to 
M. 1 := X(t, -)A4 at the point X(t, xq). If these limits exist from both sides of M. at Xq, using the Holder 
continuity of F, then the Rankine-Hugoniot condition (|1.26[) holds at the point X(t, xq). 

To prove the regularity assertions in Theorem 11.61 (a), we consider the following two case. 

Case 1. If po(xo+) = 0, then using a similar argument as that in the proof of Thcorcm ll.41 we have 
p(X(t, xo)+, t) = for each t > 0. Thus, the maps 1 1— > p(X(t, xq)+, t) and 1 1— > P(p(X(t, Xo)+, t)) are in 
C'([0 ; oc)). 

Case 2. If po(xo+) > 0, then there exist p > and r > 0, such that po\E + nB ro (x ) ^ P- Denote 
A s := A(p s ) and A := A(p). From ([Oljl . we have 

Af + u 6 • VA 5 = -F 5 - P(p 4 ) + P(p). 

We will choose a sequence of smooth test functions {(f) S ' h }s,h>o satisfying the equation <j) 5 t ' h +div((f> 5 ' h u s ) — 
0, so that 

b s ' h A 5 dx = - [ [ 4> 5,h (F s + P(p 5 ) - P{p))dxds. (5.6) 
o Jo J 

To construct (f> s ' h , we let {xh}h>o be a sequence in E + and {rh} a sequence of positive numbers such 
that Xh — > a;o and rh — > as h —> 0, and B2r h {xh) C £7+ n S ro (xo). Then, define (/A' 1 the solution of the 
equation 

b 5 t ' h + div{0 s ' h u s ) = 0, 
^ h \t=o = ^, 

where <Pq is a smooth function with support in B rh (xh), J <\>^dx = 1 and < (ffi < C h . It follows that 
(j) 5 ' h has support in X s (t, -)B rh (x h ), J (j> S ' h (x,t)dx = 1 for t > 0, and < </> S ' h < C h {T) for < t < T. 
This last assertion is a consequence of (|2.52|) . 
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In the following, we will take limits as 8 and h go to zero in (|5 . 6|) . Form Lemma 15 . 21 we have X s 
converges uniformly to X on [0,t] x B rQ (xa) as 8 — > 0. Then, we can easily prove that, for each h > 0, 
there is a <$o(/i) > such that 

X s (s,-)B rh (x h ) d X(s,-)B 2rh (x h ), 0<S<S (h). (5.7) 

We now obtain uniform bounds of the three terms in (|5.6p . For fixed t > 0, we get 

b 5 ' h A s dx = J /'''(A 5 - A)dx + J (f) s > h kdx := I + II. 

From the convergence of p s and (fn5|) , we obtain that / — > as 8 — ► 0, and J is bounded by C h (T). 
Also, from (|5.7p , we have 



essinfA(-,t)|x(t,.)B 2r . h (x h ) < ^ < esssupA(-,i)| x(t) .) B2rh ( X(i) , 
for sufficiently small 8. Thus, there exists 8o(h) > 



essmfA(-,t)\x(t,-)B 2rh {x h ) -h < J <p d ' h A d dx < esssu P A(-, t)|x(t,.)&r h (**) + fe, (5.8) 

for <5 < 5 (/i)- 

Let A(4+,t) := A(p(x t Q +,t)). Then, 

A(xp+,i) = lim esssupA(-,t)| B ,( x *)ns* = 1™ essinfA(-,t)| B (l i )n i;t . 

For each r' > 0, we can choose h r < > such that X(t, -)B 2 r h {xh) C Br^x^) n for all /i < /v- Thus, 
for each r' > and all h < h r ', 

essmfA(-,t)| Bp , (x t )n ^ < essm{A(-,t)\ X ( t ,-)B 2rh (x h ) < esssupA(-,t)\ Br , {x t o)nE t + . 

Taking first the liminf and limsup as h — * and then the limit as r' — ► 0, we obtain that 

lire i essinfA(-,t)| x(ti . )s {xh] = A(x t +,t). 

n — >u 



Similarly, we have 
From (15.81). we have 



limesssupA(-,£)|x(t,.)B 2r , ( Xh ) = A(x' +,t). 

h — *0 n 



lim lim / (<j) S ' h A & )(x,t)dx = A(x t Q +,t). (5.9) 



Similarly, we obtain 



lim lim / ((f> S ' h A 5 )(x,0)dx = A(x o +,0), (5.10) 

lim lim / [(<p S ' h F s )(x lS )dxds= [ F(x s Q +,s)ds, (5.11) 

h->0 6->0j J J Q 

lim lim / f ^ h {P(p 5 ) - P(p)}dxda = [ {P(p(x» +, a)) - P(/5)}d S . (5.12) 



From ([576]) and (|5~9| - (|5~T2)) . we have 



A(4+, t) - A(x +, 0) = - / {F(a;g+, a) + P(p(x s +, a)) - P(p)}ds. (5.13) 



Since F is locally Holder continuous in R 2 x (0,oo), standard ODE theory implies that the map t 
A(xl+,t) is in C([0,oo)) n ^((O, oo)). From (fTTT3|) . JTTH|) . (j2~lS) and (|23T|) . we have 

||F(-,s)|| L oods < C / (a- 2 -i)w<CtJW, 
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Thus, map t i-> A(x t Q +,t) is in ([0, 00)) H (^((O, 00)). From (fTTTB"]) and Theorem Ol the maps 

i 1—* p(X(t, x )+, t) and i > P(p(X(t, Xq)+, t)) are in 

Since F is locally Holder continuous in M. 2 x (0, 00), then the map t 1— > divtt(X(t, xo)+, is locally 
Holder continuous on (0, 00). These finish the proof of Theorem 11.61 (a). 

Then, we will prove Theorem 1 1.61 (c) as follows. From (|5. 13|) . we have 

J 



[A(4,t)]-[A(x ,0)] = - / [P(p(x s +,s))]ds = - a(s,x )[A(x s +,s)]ds. 

Jq Jo 

Using Gronwall's inequality, we can finish the proof of Theorem 11.61 □ 

6 Nonphysical solution 

Proof of Theorem 11.71 

First, from (|1.5[) . we have 

d r (d r v + -v) = 0, mV'(U), 
r 

where U = {(r,t) e M+ x R+\r e (a(t),b(t))}. From the regularity (fl~T5)) , we have 

v (r,t) = a (i)r + — , in X>'(t7), (6.1) 
r 

where a, f3 € ^; oe ([0, oo)). 

Then, from (|1 .30|) and the regularity of the solution in Theorem ll.il we have 



dt 

where 



E(t) = {vr ((A + 2p){v r + - P( e )) } 



r=o(t)— 

E(t)= j [^gv 2 + G(g)]rdr+ [ \\ + 2/i)(« r + y/,-. 
From (TlTTSj) and |[T^B]1 . we get 



11 \2 / Jo Jo 



From (|6.1[) , we have 



M = 0, [(A(g) + 2fi)(v r + -)- P(g)] = 0, at (a(t), t), 
and since g{a(t) + 0, i) ~ 0, we obtain 

—E(t) = a(t)v{a(t),t)(X(0) + 2p){v r + -)(a(t) + 0, t). (6.2) 
at r 

{Vr + _ )(o(t) + o, t) = 2 Q2(t) „ fe2W ■ (6-3) 

From (|6.2p - (|6.3p . we can immediately obtain (|1.33p . □ 



7 Non-global Existence of Regular Solutions 

Proof of Theorem 11.81 

Since the support of the initial spherical symmetric density po is compact, we can assume that 

supppo = 0(0) M 2 ||x| < i? }, 

for some i?o > 0. We set 

n{t) = {x = X(t,a)\a e O(0)}, t e [0,T], 
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where X(t,a) is defined in Theorem 11.21 From the transport equation (JTTTJi, one can easily show that 
suppp(x,f) = fi(i) = R 2 \\x\ < R(t)}. From fO]) . we have 

v(r,t) = a(t)r + ^-, inV'CU), 
r 

where u(x,t) = v(r,t)f, \x\ =r,U = {x £ Vl c (t),t g [0,T]}. Since u g C 1 ([0, T];H k ), wehavew(x,t) = 
in x g O c (i), and fi(t) = fi(0) for all < t < T. 

Now, we introduce the following functional as in (TTJ HB] : 

H(t) = J(x- {l + t)u) 2 pdx+ —?—^{l + t) 2 J Ap~<dx (7.1) 

x 2 pdx- 2(1 + t) [ pu ■ xdx + (1 + t) 2 [ ( pu 2 + -^-pA dx, t g [0,T]. 



7-1 

Using Q1.29[) - (|1.30|) . the Cauchy-Schwarz inequality and Holder's inequality, we obtain 
H ,( t j = 4(2- 7 )(l + t) f Ag^rdr + 4(1 + t) f \{v r + -)rdr 



7-1 

f°° v 
-2(1 + 1) 2 / (A + 2p)(v r + -frdr 



r 



< 4 ( 2 -7)(l + j~ AfPrdr + 2 j~ C( f rdr 

< 4 ( 2 -7)(l + *) f Ag J rdr + ^ f + O(0), (7.2) 

7 - 1 7o 7 Jo 7 

where i g [0, T]. From (JT7TJ)— ^T^J), we get 
and 

< (i + tf- 2 ^- 2 -^ Uo) + 2c(7 - ^ m J\l + ^-S^^A, 

< (l + t)^*r S(0) +22(7^^3^^ 

7 

where i g [0, T] and 

p W = ; ^(l-(l + 3 - 27 ), 27-3^0, 
W \ (l+t) 4 - 2 ^ln(l + i), 2 7 -3 = 0. 

From (|7.ip . we have 

^ < 1^1(1 + i) 2 " 2 ^(0) + c(7 - 1)(7 - /3) n(0) e 3sg ^F(t)(l + r 2 := G(t). 
2 7 

By the conservation of mass and Holder's inequality, we obtain 

K dx = ! pdx<( [ p>dx) ' (S!(0))^ < (^Sl) (fifO)) 3 * 1 , ( £ [0,T]. 



Since lim^oo G(t) = 0, the above inequality implies that T must be finite. □ 
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